Can decaying modes save void models for acceleration? 
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The unexpected dimness of Type la supernovae (SNe), apparently due to accelerated expansion 
driven by some form of dark energy or modified gravity, has led to attempts to explain the ob- 
servations using only general relativity with baryonic and cold dark matter, but by dropping the 
standard assumption of homogeneity on Hubble scales. In particular, the SN data can be explained 
if we live near the centre of a Hubble-scale void. However, such void models have been shown to be 
inconsistent with various observations, assuming the void consists of a pure growing mode. Here it is 
shown that models with significant decaying mode contribution today can be ruled out on the basis 
of the expected cosmic microwave background spectral distortion. This essentially closes one of the 
very few remaining loopholes in attempts to rule out void models, and strengthens the evidence for 
Hubble-scale homogeneity. 

PACS numbers: 95.36.+X, 98.65.Dx, 98.80.Es, 98.80.Jk 



I. INTRODUCTION 

In recent years there has been considerable interest in 
using observations to try to confirm the assumption that 
the Universe is homogeneous and isotropic on the largest 
scales. While isotropy about our own worldlinc follows 
almost directly from the observed isotropy of the cos- 
mic microwave background (CMB), radial homogeneity 
is difficult to confirm on Hubble scales. Galaxy surveys 
have found that the distribution of luminous matter to 
redshifts of at least z ~ 0.5 is largely consistent with 
the standard cosmological constant plus cold dark matter 
(ACDM) model (e.g., see [1]). However, radial inhomo- 
geneity on the largest scales could be difficult to disen- 
tangle from rcdshift-dcpcndcnt effects such as evolution, 
so it is not clear how much current surveys actually tell 
us about these scales. 

Given the isotropy of the CMB, the possibility of ra- 
dial inhomogeneity might appear very contrived: it would 
imply that we are near the centre of a (nearly) spheri- 
cally symmetric Universe, hence apparently violating the 
Copcrnican principle. However, just such a situation has 
received much attention recently, in the context of late- 
time acceleration. It was realized some years ago [2— 
4] that if we are situated in an extremely large spher- 
ically symmetric undcrdensity, or void, the luminosity 
distance-redshift relation of Type la supernovae could be 
explained without the need for a cosmological constant 
or dark energy, or a modification of gravity. (An ear- 
lier study of the potential pitfalls of living in a spheri- 
cal undcrdensity while assuming a homogeneous model 
can be found in [5].) These models exchange the coinci- 
dence problem of ACDM for a violation of the Copernican 
principle (although they are not free of temporal tuning 
problems [6], nor do they address the "old" cosmologi- 
cal constant problem). For a recent review of these void 
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models and the various constraints on them, see, e.g., [7]. 

General geometrical tests of large-scale homogeneity 
have been proposed (e.g., see [8, 9]). In the context of 
void models for acceleration, the strongest and most ro- 
bust results come from the requirement of fitting the full 
power spectrum of CMB temperature anisotropics. It has 
been shown that growing mode void models which do fit 
the CMB predict a local Hubble rate far lower than obser- 
vations indicate [10-14]. Claims have been made [15] that 
early radiation inhomogeneity can provide a loophole to 
this conclusion, although this appears unlikely due to the 
free streaming and rapid redshifting of the radiation [11]. 
Modifying the primordial perturbation spectrum might 
provide a loophole [10, 16], although order- unity depar- 
tures from scale invariance, and substantial fine tuning, 
would be required [11]. 

Particularly promising tests of homogeneity are those 
that rely on the scattering of light from inside our past 
light cone. Void models generically predict large CMB 
dipoles for off-centre comoving observers [5], and hence 
should generate anisotropics via the kinetic Sunyaev- 
Zel'dovich (kSZ) effect [17] (as first suggested in [18]). 
In Rcfs. [19, 20], this effect was used to put constraints 
on void models using galaxy cluster observations, and 
in [21] the linear kSZ effect due to all structure was 
found to essentially rule out a large class of void models. 
In [22], it was shown that strong linear kSZ constraints 
on (growing-mode) void models persisted under a fully 
relativistic treatment of the problem using the Lemaitre- 
Tolman-Bondi (LTB) spacetime [23-25]. However, it was 
stressed in [22] that there were many caveats to the lin- 
ear kSZ calculations, both technical and relating to the 
poor state of model-independent knowledge of the local 
baryon fraction and free electron density power. 

Another method, which is much less susceptible to 
the uncertainties which plague the kSZ approach, in- 
volves observations of the CMB spectral distortions due to 
Compton scattering from inside the light cone [18]. While 
the kSZ approach relies on the knowledge of the free elec- 
tron perturbation power, the calculation of the Compton 
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y distortion requires only background information. Al- 
though tight constraints on growing-mode void models 
using the Compton y distortion have been presented [26] , 
more recent studies using a fully consistent LTB approach 
have found only very weak constraints [11, 22]. 

Even though the kSZ effect and spectral distortions 
both probe the inside of our past light cone, they have 
so far only been used to gauge possible inhomogeneity 
on the surface of the light cone, at relatively late times. 
This is because of the usual restriction to growing modes 
of the LTB solution. Both growing and decaying modes 
arc possible [27], but decaying modes are usually dis- 
carded on the basis that they imply extreme inhomo- 
geneity at early times, and hence appear difficult to rec- 
oncile with standard inflationary scenarios [28]. The as- 
sumption of vanishing decaying modes also enables tests 
based on structure [10, 11, 28], since the growing-mode 
void can itself be treated as a linear perturbation from a 
Friedmann-Lemaitre-Robertson- Walker (FLRW) space- 
time at early times. The absence of decaying modes 
means that the large dipolcs in the LTB models usually 
studied are mainly local to the scatterers — they can be 
considered as the result of "peculiar velocities" with re- 
spect to the CMB frame — and hence the ability of these 
techniques to probe the inside of the light cone has not 
been fully exploited. 

In standard cosmological models, any decaying modes 
are assumed to have negligible amplitudes by the end 
of inflation. Decaying modes on our last scattering sur- 
face (LSS) have been obscrvationally constrained to be 
a subdominant component [29] (see also Rcf. [30] for an 
estimate of the effect of small-amplitude decaying modes 
on the LSS). In the context of void models for accelera- 
tion, however, a decaying mode may be localized around 
the observer (near the centre of spherical symmetry) and 
hence not extend to the LSS. Nevertheless, such a de- 
caying mode would be visible to a scatterer at the ap- 
propriate rcdshift, and so should be expected to generate 
substantial spectral distortions and kSZ anisotropics. 

It would be ideal to rule out decaying modes in void 
models on the basis of observations, rather than theoret- 
ical prejudices. In addition, it has been suggested that 
the extra freedom from incorporating a significant de- 
caying mode contribution today might ease the present 
severe constraints on void models [30, 31]. Therefore, 
in this work, full use of the power of the y distortion 
to probe inside the light cone is made to study the ob- 
servational implications of decaying modes and hence to 
constrain inhomogeneity at the earliest times. The y dis- 
tortion is studied rather than the kSZ effect because of 
the above-mentioned ambiguities with the kSZ approach. 
The simpler y distortion calculations, combined with the 
current observational upper limits, will, nevertheless, be 
sufficient to rule out any contribution of decaying modes 
significant today, for all but extremely narrow decaying 
mode profiles. 

In Sec. II, the LTB model is introduced and vari- 
ous properties of decaying modes discussed. Next, in 



Sec. Ill, calculation procedures for the y distortion are 
introduced, both at nonlinear order and under the linear 
approximation. Results are presented in Sec. IV, before 
conclusions are made in Sec. V. A covariant derivation 
of the LTB null geodesic equations is presented in the 
Appendix. Units are chosen such that c = 1. 



II. LTB DECAYING MODES 
A. General LTB solution 

As mentioned in the Introduction, the general LTB 
spacetime consists of both a growing and a decaying 
mode [27]. This section begins with a review of the iso- 
lation of these modes (based on Ref. [28]), followed by a 
description of some technical properties of the decaying 
modes. 

For a spherically symmetric distribution of pressure- 
less matter, Einstein's equations can be solved exactly, 
resulting in the LTB solution. The metric can be written 



ds 2 



-dt z 



Y' 



1 - K 



-dr + Y dCl , 



(1) 



where a prime denotes the derivative with respect to co- 
moving radial coordinate r, and t is the proper time along 
the comoving worldlines. The function K = K(r) < 1 
is arbitrary, and the areal radius Y = Y(t,r) is given 
paramctrically by 
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< K < 1. 



Here M(r) is a free radial function, which is set to 



M(r) 



as a gauge condition (this implies that we 



cannot extend the solution past the "equator" of a closed 
LTB model). 

Various physical quantities can be expressed in terms 
of these functions (see, e.g., [28]). In particular, 
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where p, (9, and £ arc the matter density, expansion, and 
shear scalar, respectively, for the comoving worldlines. 
For comoving shear tensor cr M „ and radial spatial unit 
vector r M (r^u^ = 0), we have £ = a^r^r" . is 
the Ricci curvature of the spatial comoving-orthogonal 
hypersurfaces. The radial and transverse expansion rates 
are given, respectively, by 

Y' Y 
ff|| = y7= H ± = y , (8) 

where the overdot denotes the derivative with respect to 
t. 

This exact solution contains another free radial func- 
tion, t B = t B (r), which is known as the "bang time" 
function since t = t B (r) implies Y = 0, leading to diver- 
gences in each of p, 8, and ^R. In the models of interest 
here, this will correspond to the cosmological singularity. 
When t B (r) is not constant, it is often stated that the 
big bang occurs at different times at different locations. 
However, we are of course free to define a new time coor- 
dinate t such that the big bang occurs "simultaneously" 
at t = 0. The physical content of a varying t B is related 
to the divergence of covariant quantities such as the shear 
scalar, £, as t — > t B . 

As described in [28] , small-amplitude variations in spa- 
tial curvature, ^R, correspond to growing- mode pertur- 
bations about an Einstein-de Sitter (EdS) FLRW model, 
while appropriately small variations in bang time cor- 
respond to decaying modes. Therefore, according to 
Eq. (7), to isolate a pure decaying mode we must set 
K(r) = 0. The LTB solution, Eqs. (2) and (3), then 
becomes 

r-(f)V«.r. W 

In the remainder of this work, only pure decaying modes 
will be considered. Justification for this will be discussed 
in Sec. V. The main motivation of this study is to ex- 
amine the effect of decaying modes that are significant 
today. In Sec. IIIC we will see that this implies that 
we are interested in bang time fluctuations St B (r) ~ to, 
where to is the proper age today. 



B. Properties of decaying modes 

Since significant decaying modes are not normally 
expected in standard inflationary scenarios, it will be 
worthwhile to explicate some of their counterintuitive 
properties. Indeed, in the years before inflation was pro- 
posed, decaying modes were studied extensively, as there 
was no reason to ignore them a priori. In particular, un- 
der the guise of "delayed cores" of the big bang, they were 
proposed as an explanation for quasars [32]. LTB de- 
caying modes are generalizations of Schwarzschild white 
holes, which are time-reversed analogs of black holes. 



An important property of LTB models involves the 
presence of shell crossings, where comoving worldlines in- 
tersect and hence the pressure-free assumption becomes 
invalid. In [33], criteria were derived to ensure that no 
shell crossings occur; for the case of our K(r) = pure 
decaying mode spacetime, they become 

t' B < 0. (10) 

This condition restricts us to situations in which the big 
bang singularity is "delayed" at the origin, i.e. t B (r = 
0) > t B (r). Throughout this work, the condition (10) 
will be assumed to hold. 

Another relevant feature of LTB decaying modes is the 
generic presence of cosmological blueshiits, and possibly 
divergent blueshifts. In [34] it was shown that light rays 
emitted near the initial singularity will generically be 
blueshifted at late times when t' B ^ 0. This might mean, 
e.g., that an observer whose LSS intersected a decaying 
mode would observe the CMB temperature in the direc- 
tion of the decaying mode to be greater than the actual 
recombination temperature! In addition, in [35] it was 
pointed out that under certain circumstances, divergent 
blueshifts can arise in such models, with the result that 
the white hole will be unstable and will convert into a 
black hole. This blueshift behaviour can be understood 
broadly by noting that for t' B ^ 0, the spacetime be- 
comes shear dominated at early times. Since the shear is 
defined to be trace free, this means that in some direc- 
tions the comoving worldlines will be contracting near the 
singularity. In particular, the radial expansion rate H\\ 
will be negative close enough to the singularity, leading 
to blueshifts for null rays that are approximately radial. 
(A related discussion appears in Ref. [36].) 

In the void model context, a decaying mode near the 
centre would not be directly observable by a central ob- 
server. However, a scatterer at the appropriate red- 
shift down the central observer's light cone would see the 
blueshifts due to the decaying mode, and hence would 
observe a strongly anisotropic CMB sky. This suggests 
that the scatterer would produce significant CMB spec- 
tral distortions (or kSZ anisotropies). Of course in a re- 
alistic cosmology, the dust source approximation of LTB 
breaks down at early times, when radiation becomes im- 
portant. It is not clear how the blueshift behaviour will 
be modified in the radiation era. However, since recombi- 
nation occurs somewhat after radiation domination, the 
large blueshifts will be incurred during the dust era, when 
the LTB solution is reliable. Thus the constraints from 
the spectral distortion will be sound. 

III. CALCULATING THE SPECTRAL y 
DISTORTION 

A. The y distortion 

The y distortion arises from the Compton scattering 
of anisotropic CMB radiation from inside our past light 
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cone into our line of sight. In the single-scattering ap- 
proximation, and in a spherically symmetric spacetime, 
it can be written as [37] 



Ye / ° dZs Jz I d £ sin £(i + cos2 ln 



21 t(z s ,o 



T{z s 



(11) 

Here r is the optical depth, and the outer integral only 
extends to the redshift of (the assumed abrupt) reion- 
ization, z TC . T(z s ,£) is the CMB temperature seen by a 
scatterer at redshift z s from the central observer. Spher- 
ical symmetry implies that the temperature is a function 
of only one angle, £. A radially outwards directed ray 
is chosen to have £ = 0. Although in standard mod- 
els we have z re — 10, the reionization redshift could be 
somewhat different in void models. Nevertheless, in the 
models studied here, most of the contribution to the red- 
shift integral will come from z s ~ 3, so uncertainty in z rc 
should not affect the final results significantly. The factor 
d,T/dz s can be readily written in terms of the background 
quantities as [22] 



dr _ <TT fb (2 - Yhc) Pm{z s ) 
dz s 2m p (l + z s )H\\ (z a ) 



(12) 



Here or is the Thomson cross section, ft, = pj,/ Pm is the 
baryonic to total matter fraction, Ijje is the helium mass 
fraction, and m p is the proton mass. 

In the case that the fluctuations over the sky in T(z s , £) 
are small, and dominated by the dipole /3(z s ), Eq. (11) 
becomes [11] 



y= wJ dZs dz-/ (Zs)2 - 



(13) 



We will see that this is generally not a good approxima- 
tion for decaying mode void models. 

The best current constraints on the y distortion come 
from the COsmic Background Explorer (COBE) satellite, 
which found y < 1.5 x 10~ 5 at 2a confidence [38]. 



B. Nonlinear calculation of CMB anisotropics 

In order to evaluate the y distortion expression, 
Eq. (11), we need to know the CMB temperature 
anisotropy seen by a scatterer at redshift z s , T(z s , £). To 
evaluate this temperature, we need to determine the red- 
shifts along various (generally nonradial) null geodesies. 
The general exact expression for redshift (e.g., see [39]), 
applied to a geodesic comoving congruence of worldlines, 
gives 



1 + z(t,£) = exp 



+ <r fl „n fl n l/ dt 



(14) 



Here z(t, £) is the total redshift incurred along the null 
geodesic between an arbitrary point, at proper time t, 



and the ray's endpoint, at t = t e . The integral is eval- 
uated along the null ray. The vector n^(t) is the nor- 
malized projection of the geodesic's tangent vector, v 1 *, 
orthogonal to the comoving timelike vector field it M (see 
the Appendix). In other words, n^(t) is the spatial di- 
rection of propagation of the null ray seen by a comoving 
observer. At the endpoint of the geodesic, n M (t e ) is di- 
rected at angle £ with respect to the radially outwards 
direction, i.e. r^n^(t e ) = cos£. 

In order to evaluate Eq. (14), we must first know how 
to determine the past-directed null geodesic with initial 
condition (IC) £ or n^(t e ), and then we must know where 
to stop evolving the geodesic, i.e. we must know how to 
find the LSS. The first problem is straightforward, in that 
we must solve the null geodesic equation 



v?„v v = 0. 



(15) 



Extracting various components of Eq. (15) (see the 
Appendix) gives 



(16) 
(17) 
(18) 
(19) 



dt 




dX 


= 7, 


dr 
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dX 


~ Y > ' 


d6 
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dX 


~ y2' 


dn r 


1 , 2 


~dX 


= 27K 



Here 9 is the standard spherical angular (comoving) co- 
ordinate of the photon (the second spherical angle can 
be chosen to be constant), A is an affine parameter along 
the null ray, 7 is defined by Eq. (16), n r = n^r^, and L 
is a constant. Finally, Eq. (14) now gives 



*(*,£) = -■ 

7e 



(20) 



The set of coupled ordinary differential equations (16) 
to (19), together with the ICs (i.e. the endpoint values) 
t e , r e , 9 e = 0, n r , e = cos£, and j e , can then be solved 
numerically with standard techniques. 

Finding the LSS, i.e. knowing where to stop the 
geodesic integration, necessarily involves some assump- 
tions, since the details of recombination in the region of 
the decaying mode are unclear. In this work, it will be 
assumed that recombination occurs at the same energy 
density and local temperature within the decaying mode 
as asymptotically outside. Schematically, the approach 
is to first use Eqs. (16) to (20) to evolve a null geodesic 
from the centre of symmetry today into the past to a 
redshift of zls = 1091, where the energy density pls is 
evaluated. Then, for each scatterer redshift z s , several 
null rays are propagated, again into the past, from the 
scatterer in several spatial directions as characterized by 
the initial angle £. These null rays are evolved until the 
local energy density reaches /?lSj at time £ls(£)i which 
is taken to define the LSS in that direction. The total 



5 



rcdshift, z(iLS (£),£)' incurred from the central point to- 
day, to the scatterer, to the LSS, is then evaluated using 
Eq. (20). Finally, the temperature observed at the scat- 
terer in direction £ can be written 



T(z s ,0=T 



(l + z L s)(l + z a ) 

1 + *(*LS(0,£) 



(21) 



where To is the (unscattered) CMB temperature today. 
This expression can then be used to evaluate the inte- 
gral for the y distortion, Eq. (11). Note that this tech- 
nique generalizes that of Ref. [11] to the case of nonra- 
dial geodesies, and to a LSS determined by fixed density, 
rather than fixed proper time (which is a reasonable ap- 
proximation in the case of vanishing decaying modes). 
The dipole /3(z s ) in Eq. (13) is calculated in the LTB 
model by propagating past-directed null rays radially in- 
wards and outwards from the scatterer to plSj and com- 
paring the resulting redshifts as described in [11]. 

Importantly, the uncertainties about recombination in 
the decaying mode region will likely affect the calcula- 
tions of T(z a ,£) and hence of y. Since the temperature 
at recombination is determined by local atomic physics, 
the assumption that the recombination temperature is 
constant is likely reasonable. However, it is certainly pos- 
sible that there may be order-unity uncertainties in the 
density at recombination, due to the presence of shear 
close to the singularity. (Likewise, as mentioned earlier, 
the radiation component absent from the LTB solution 
is expected to affect the T(z s ,£) calculation at the level 
of tens of percent.) Nevertheless, we will see that the 
results presented here are powerful enough that order- 
unity variations in the density will not affect the final 
conclusions. 



C. Linear calculation of CMB anisotropics 

As a check of the nonlinear calculation described in 
Sec. IIIB, and especially given that such calculations for 
decaying modes have apparently not been made previ- 
ously, it will be useful to perform an alternative calcu- 
lation of T(z s ,£) in the case where a linear description 
is valid. A decaying mode associated with bang time 
fluctuation Stsir) will be accurately described by lin- 
ear theory at time t when StB{r)/t -C 1 [28]. (This ex- 
plains the condition Stsir) ~ to for decaying modes to 
have significant amplitude today.) Therefore, if we satisfy 
^bMAls *C 1 at the time of recombination, then the 
spacetime will be accurately described by a linear per- 
turbation from FLRW for all times after recombination 
(since growing modes are ignored). 

In particular, we can describe a small-amplitude decay- 
ing mode in terms of linear metric perturbations in some 
gauge. A convenient choice is Newtonian (or zero shear) 
gauge, for which the vanishing of anisotropic stress im- 
plies that there is only one unique metric perturbation, 
ip, which describes both the lapse and isotropic spatial 
metric (curvature) perturbations (e.g., see [40]). Then 



the temperature anisotropy viewed in direction is de- 
scribed by the Sachs- Wolfe (SW) effect [41]. In the ap- 
proximation of abrupt recombination (which is a good 
approximation on large scales) , this can be written [39] 



4'dt + -ip 



'tLS 

2 1 
3 



9 \H 







i 2 H 2 



U + Hip 



(22) 



Here a and H are the FLRW background scale factor and 
Hubble rate, and the quantities outside the integral are to 
be evaluated at the point on the LSS in direction . T e 
is the background temperature at the null ray endpoint 
(normally taken to be the scatterer at z s ). (Note that 
Eq. (22) corrects a typographical sign error in [39].) 

Equation (22) is general in that it applies to growing 
or decaying modes. However, we can simplify it for the 
case of decaying modes. In a dust background, the metric 
perturbation ip satisfies the equation (e.g., see [40]) 



ip + AH'ip = 0. 
The decaying mode solution is 



5/3 



which implies 







-Hip. 



(23) 



(24) 



(25) 



Substituting Eq. (25) into Eq. (22) gives 



6T(n>- 



-5 / ip(t,r(t))H(t)dt + 2ip 

2 V 2 1 

— s — o ip H — iu/. 

9a 2 H 2 ^ H 



(26) 



In an EdS background, the radial component of the null 
trajectory is given by 



r(t) = \/r 2 + Ar 2 — 2r e Arcos£, 



where 



Ar(t) 



a e H e 



1/3- 



(27) 



(28) 



Equation (26) is the final expression for the decaying 
mode SW effect. Note, in particular, that it includes 
an integrated SW component, since even in an EdS back- 
ground, ip is time dependent. 

The final step is to relate the metric perturbation ip(r) 
to the bang time fluctuation, t B {r). To do this, first note 
that the relation 



Sp _ 2 V 2 
~~p ~ 3a 2 H 2 



1>, 



(29) 
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where Sp is the comoving density perturbation, holds 
for both growing and decaying modes. Then, expand- 
ing Eq. (4) with Eq. (9) in terms of the small parameter 
tB(r)/t, we can write 



Sp 
P 



t_B_ 
t 



t' B M 



t M 



TT7 +0 



(30) C 



Note that this result generalizes the corresponding result 
in [28] to the case of significant bang time gradients, t' B . 
Now, given a bang time function, all the pieces are in 
place to perform the linearized calculation of the CMB 
anisotropics, and hence of the y distortion. 



IV. RESULTS 
A. Bang time profile 

In order to calculate the y distortion, we must first 
choose a bang time function. A convenient choice is a 
Gaussian profile, specified by 



*B(r) = tB, m e 



-r 2 /L 2 



(31) 



with amplitude ts.m > and characteristic width L. 
This specification is all we need to calculate the y distor- 
tion according to the nonlinear prescription of Sec. IIIB. 
However, for the linear calculation of Sec. IIIC, we first 
need to determine the corresponding metric perturba- 
tion, ip. It is straightforward to show that Eqs. (29) and 
(30) arc satisfied when ip takes the form 



(a LH y 



tB- 



2*0 



5/3 



-r 2 /L 2 



(32) 



Here the right-hand side has been written in terms of the 
quantities today, do, Hq, and *q, since they are constant. 
The fact that ij>(t,r) is also a spatial Gaussian explains 
why the Gaussian choice for *s(r) is convenient. 

For the calculations presented in this section, units 
were chosen such that 100 km s -1 Mpc -1 = 1 and AttG = 
1. This means that realistic values of the Hubble pa- 
rameter today are H ~ 0.7, and that values of the 
matter density and proper time today are both of order 
unity. It is difficult in general inhomogeneous cosmolo- 
gies to unambiguously describe distances. Therefore, the 
bang time profile width L will be translated into a corre- 
sponding redshift down the central observer's light cone, 
zl = z(r = L). As already mentioned, we will be pri- 
marily interested in the regime tB,m ~ to, in which the 
decaying modes will have a significant effect at late times. 
All of the calculations in this study used a value H =0.7 
and integrated to a reionization redshift of z Te = 11. 



B. Convergence to linear theory 

Before calculating the y distortion, it will be important 
to check that the nonlinear calculation of temperature 
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FIG. 1: CMB temperature anisotropy visible to a scatterer at 
(f , r) = (1.0, 1.15). The angle £ is measured from the radial di- 
rection, i.e. £ = corresponds to the centre of the anisotropy. 
Each solid curve represents a nonlinear calculation using the 
method of Sec. IIIB, while each dashed curve used the linear 
method of Sec. IIIC. The bang time amplitude took the val- 
ues t B , m = 1-0 x 10~ 3 , 1.0 x 10~ 4 , and 1.0 x 10~ 5 , for the 
top (blue) curves, middle (red) curves, and bottom (black) 
curves, respectively. The nonlinear and linear calculations 
converge at the smallest bang time amplitudes, as expected 
(the bottom two curves are almost indistinguishable). 



anisotropics agrees with the linear calculation, in the ap- 
propriate limit. In Fig. 1, the anisotropy T(f s ,r s ,£) is 
displayed for a scatterer at position (t s ,r s ) = (1.0, 1.15), 
calculated using both the nonlinear and linear ap- 
proaches. The scatterer's position is such that the scat- 
terer's past light cone intersects the decaying mode some- 
what away from the peak of the Gaussian (although the 
results are qualitatively similar for any scatterer posi- 
tion). The calculations used a profile width of L = 0.25, 
which corresponds to zl — 0.35. As the decaying mode 
amplitude decreases, it is apparent that the two 

calculations converge, as expected. This is particularly 
reassuring, considering how very different the anisotropy 
calculations are for linearized EdS and nonlinear LTB. 
In addition, notice that the anisotropy becomes of order 
unity already for the relatively small amplitude ts.m = 
1.0 x 10 -3 . For decaying modes of significant amplitude 
today, i.e. tB,m ~ *o ~ lj we therefore expect extreme 
anisotropies at the scatterers, due to the blueshifts men- 
tioned previously. 

For the calculations in Fig. 1, the time of last scat- 
tering corresponds to *ls = 2.8 x 10 -5 (in the linear 
regime, in which £ls is independent of position). There- 
fore, the three pairs of curves in that figure correspond 
to tB, m /t~LS = 36.1, 3.61, and 0.361, from top to bottom. 
We can see that the linear approximation is good even 
for ratios £_B,m/iLS of order unity. Thus the condition 
tB.m/ths <§C 1 given in Sec. IIIC for the validity of linear 
theory is too stringent. To understand the reason for this, 
note that the anisotropies shown in Fig. 1 are small, even 
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when tB,m/thS is of order unity. This occurs because of 
a near cancellation of the first two terms in the decaying 
mode SW effect, Eq. (26). Therefore, a metric pertur- 
bation ip of order unity can produce small anisotropy, 
ST/T. This is in contrast to the familiar growing mode 
case, where the SW anisotropy will be of the same order 
as if). Intuitively, a decaying mode decays, and hence has 
much less effect integrated over the null ray than a grow- 
ing mode. Note, finally, that this also means that it is 
important to include the gradient terms in Eq. (26), even 
for the Hubble-scale decaying modes considered here. 



C. y distortion 

With all of the tools in place, we can now calculate the 
y distortion. Figure 2 displays the y distortion as a func- 
tion of the dimensionless bang time profile amplitude, 
tB.m/to- These calculations also used a profile width of 
L = 0.25, which corresponds to zl ~ 0.35. Curves are 
shown for the full-sky integration based on Eq. (11), to- 
gether with the dipole approximation, Eq. (13). Also, 
anisotropics are calculated using the nonlinear approach 
of Sec. IIIB, as well as the linearized SW approach from 
Sec. Ill C. It is evident, again, that the nonlinear calcu- 
lation approaches the linear one for small bang time am- 
plitudes, and that both exhibit the expected quadratic 
dependence of y on tB, m /to in the linear regime. It is 
clear as well that the dipole approximation substantially 
overestimates the spectral distortion, especially at the 
lower amplitudes. Most striking, however, is the fact 
that the y distortion exceeds the COBE limit by orders 
of magnitude, for decaying modes which are significant 
today (i.e. for tB, m /to ~ 1), at least for decaying modes 
with width z L — 0.35. 

We should expect that the y distortion will decrease 
as the decaying mode profile width decreases, due to 
the smaller solid angle sourcing the angular integral in 
Eq. (11). Also, the anisotropics themselves should be af- 
fected [e.g. via the L 2 factor in Eq. (32)], although it is 
not clear exactly what L dependence is expected, due to 
the subtle cancellations that occur in the SW expression, 
Eq. (26). In Fig. 3, the y distortion is plotted versus 
decaying mode width, z L , at fixed bang time amplitude, 
tB,m/to = 0.2. These calculations used the nonlinear 
approach of Sec. IIIB. A marked decrease in y as z L 
decreases is visible, with y dropping below the COBE 
limit for z L ^ 0.02. Also notable is the fact that the 
dipole approximation severely underestimates the drop 
in y. This is understandable, since the dipole approxi- 
mation misses the solid angle effect mentioned above. Fi- 
nally, note that, for the profiles able to evade the COBE 
constraint (zl 0.02), the profile width must be very 
small relative to cosmological scales. 
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FIG. 2: y distortion versus dimensionless bang time ampli- 
tude, tB,m/to, at fixed profile width L = 0.25, which cor- 
responds to zl — 0.35. The black (lower) curves used the 
full-sky integration, Eq. (11), while the blue (upper) curves 
used the dipole approximation, Eq. (13). Each solid curve 
used the nonlinear anisotropy calculation of Sec. IIIB, and the 
dashed curves used the linearized SW approach of Sec. IIIC. 
The nonlinear and linear calculations converge at the smallest 
bang time amplitudes, as expected. The dotted line indicates 
the COBE upper limit of y < 1.5 x 10~ 5 at 2a confidence [38]. 



10 




FIG. 3: y distortion versus decaying mode width, zl, at fixed 
bang time amplitude, ts, m /to = 0.2. The black (solid) curve 
used the full-sky integration, Eq. (11), while the blue (dashed) 
curve used the dipole approximation, Eq. (13). The dotted 
line indicates the COBE upper limit of y < 1.5 x 10~ J at 2a 
confidence [38]. 



V. CONCLUSIONS 

As discussed in Sec. II B, cosmological decaying modes 
suffer from a variety of problems. They should have de- 
cayed to negligible amplitudes by the end of inflation, 
and no mechanism has been proposed for the genera- 
tion of decaying modes so large as to still be significant 
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today. They must satisfy the condition t' B < to avoid 
shell crossings, although there is something oddly restric- 
tive about this criterion: clearly a generic profile in the 
absence (or even in the presence) of spherical symmetry 
will violate it somewhere. Of course, without an explicit 
mechanism of formation, not much else can be said. The 
possibility of blucshift instability would also need to be 
addressed for these models. 

In addition, if we are to save void models for accelera- 
tion by adding a significant decaying mode contribution 
today, then the ICs would need to be exquisitely finely 
tuned. In the linear regime, the metric perturbation i/j 
is constant in the growing mode, but obeys Eq. (24) for 
the decaying mode. Therefore the ratio between growing 
and decaying modes increases like i 5 / 3 . A ratio of order 
unity today implies that the growing modes must be sup- 
pressed by a factor (£lsAo) 5 ^ 3 ~ 10~ 8 at l ast scattering. 
Of course this estimate will be quantitatively affected by 
the breakdown of linear theory at early times, but the 
conclusion will still hold qualitatively: In a model with 
significant decaying modes today, the extremely inhomo- 
geneous early Universe must be pure decaying mode to 
very high precision. 

One further problem with decaying mode models is re- 
lated to the shear domination at early times. Since the 
details of big bang nucleosynthesis (BBN) depend on the 
expansion history of the early Universe, we can expect 
substantial departures from standard BBN in these mod- 
els. Since the decaying mode profile must be centred near 
our own worldline, we would not expect the observed lo- 
cal light-clement abundances to agree with predictions 
based on the standard ACDM model. 

Nevertheless, in this study decaying modes have been 
given the "benefit of the doubt," in that it has been 
assumed that all of these problems could be somehow 
overcome, and the observational consequences of decay- 
ing modes have been examined. As mentioned in the 
Introduction, while growing mode LTB void models for 
acceleration have been essentially ruled out on the ba- 
sis of CMB + Hq, kSZ, and other observations, it might 
be possible that a significant contribution today of LTB 
decaying mode could allow such models to survive those 
tests. The generic presence of blucshifts near the cos- 
mological singularity of a decaying mode, however, sug- 
gests that these inhomogeneous bang time models might 
be susceptible to probes of the interior of our past light 
cone. Therefore, a procedure for calculating CMB tem- 
perature anisotropics for a scatterer whose LSS crosses a 
decaying mode was developed and used to calculate the y 
distortion in these models. Although the results will nec- 
essarily be dependent upon the uncertainties regarding 
recombination inside a decaying mode, the constraints 
were so strong that even order-unity uncertainties in the 
calculated anisotropies will have negligible effect on the 
conclusions. Decaying modes which are significant today 
and wider than Zi ~ 0.02 are ruled out by the y dis- 
tortion test. Narrower profiles might survive this test, 
but would almost certainly suffer from obvious problems 



in local surveys: the region inside zl — 0.02 would be 
much younger than the outside, and hence would have a 
substantially higher density and expansion rate. In addi- 
tion, structure would presumably have substantially less 
power on the inside. The local Universe to z ~ 0.02 is 
well surveyed, and no such discontinuities are apparent. 

It is worth stressing that decaying modes with suffi- 
ciently small amplitudes, i.e. i^m/^o 10 -5 , are not 
ruled out by the y distortion test. Although such pertur- 
bations may suffer some of the problems with decaying 
modes listed above, it is conceivable that they play a 
role in cosmology [29]. However, the effects of decaying 
modes with such small amplitudes will be negligible to- 
day, and hence will not ease the severe constraints on 
growing-mode void models for acceleration. 

Note that the philosophy of this study has been to 
examine the spectral distortions due to pure decaying 
modes. Of course one could argue that a model which 
evades the CMB + Hq and kSZ tests might consist of a 
superposition of decaying and growing modes. However, 
given that pure growing-mode models typically produce 
values of the y distortion below the COBE limit [22], it 
appears very likely that the addition of growing modes 
of typical amplitudes will not change the conclusions of 
this study significantly. Of course the combined effect 
of growing and decaying modes will generally not be a 
simple linear superposition, due to the nonlinearity of 
general relativity. But the extremely large anisotropies 
in decaying mode models are due fundamentally to the 
extreme inhomogeneity at early times, when a growing- 
mode contribution will be negligible. 

Therefore, it appears that decaying modes do not pro- 
vide a loophole to the conclusion that void models can- 
not explain acceleration. This result strengthens our 
confidence in the standard paradigm of a homogeneous, 
isotropic, accelerating cosmology, driven by some myste- 
rious form of dark energy of modified gravity. 
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Note added. — After this work was essentially complete, 
a closely related paper appeared [14], which also studied 
the observational consequences of decaying modes which 
are significant today. There are a few major differences 
to our approaches. First, Ref. [14] did not exclude the 
"advanced" big bang case, t' B > 0. This case was ex- 
cluded in the present study since the consequent early 
shell crossings were argued to render the early-time LTB 
solution, and hence the temperature anisotropy calcula- 
tions, unreliable. Another difference is that Ref. [14] used 
the kSZ effect in clusters as the basis of its constraints, as 
opposed to the y distortion used here. However, Ref. [14] 
used the dipole approximation for the anisotropies seen 
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by the scatterers. Finally, the models studied in [14] in- 
cluded LTB growing modes. The strong findings of [14] 
are in agreement with those of the present study, and 
hence reinforce the conclusion that models with signifi- 
cant decaying modes today are not viable. 



Appendix: Covariant null geodesic equations 

1. 1 + 1 + 2 formalism 

The nonradial null geodesic equations in the LTB 
spacetime have been derived many times before (e.g., 
see [34, 42-44]). In this appendix, a novel derivation 
is presented which is based on the covariant 1 + 1 + 2 
formalism, and, in particular, makes no use of Christoffcl 
symbols. 

This subsection will begin with a brief summary of the 
covariant 1 + 3 formalism. More details can be found in 
the reviews [45, 46]. This formalism is based upon a fun- 
damental timelike vector field, w M (normalized according 
to u^u^ = — 1), which is conveniently taken to be tangent 
to the dust comoving worldlines in the LTB spacetime. 
The tensor 



(A.l) 



projects orthogonal to . We can use h^ v to define a 
spatial covariant derivative according to 



D T 



h° n T 

(A.2) 

for any tensor T„ 1V3 ... Vn orthogonal to u M in all of its 
indices. 

For the case of twist-free dust, which is the case for the 
LTB model, we can decompose the covariant derivative 
of Ufj, according to 



l H;i/ 



(A.3) 



The scalar 9 measures the local volume rate of expansion 
of the congruence u M . The trace- free, symmetric tensor 
a^v measures the local rate of shear of the congruence, 
and is orthogonal to in both of its indices. For a 
dust source, the acceleration of the comoving worldlines 
vanishes, i.e. u^- p u p = 0. 

Under spherical symmetry, each comoving-orthogonal 
slice contains a preferred spacelike congruence with radial 
tangent vector r M , where r^r^ = 1. In this case, the 1 + 3 
covariant approach can be generalized to the so-called 
1 + 1 + 2 approach, which incorporates both m m and r M as 
fundamental fields [47-49]. By analogy with the tensor 
h^ v defined in Eq. (A.l), which projects into the slices, 
we can define a tensor „ by 



s v — h „ r^r v , 



(A.4) 



which projects into the two-spheres (called sheets) or- 
thogonal to both u M and r M . 



Under spherical symmetry, a spatial three-tensor such 
as the shear a^ v must be expressible in terms of r M , ', 
and a two-scalar E. Explicitly, we can write 



1 



(A.5) 



This implies that the shear scalar can be written £ 



'lit/ • 



1. Null geodesies 

The tangent vector to a null geodesic can be de- 
composed into components parallel and orthogonal to the 
dust four-velocity u p according to 



(A.6) 



where the spatial propagation direction satisfies 
n^Uft — and n^n^ — 1. By virtue of the relation 
v^Uf! = —7, the photon energy (or blackbody spectrum 
temperature) is proportional to 7. For any null geodesic 
in a spherically symmetric background, n M will be con- 
strained to a two-dimensional subspace of each comoving 
slice. In other words, n M can be decomposed into compo- 
nents parallel to r M and parallel to an angular direction, 
denoted 6^, where 6"% = 1 and 6"% = 0% = 0. Ex- 
plicitly, 



= n r r M + 



(A.7) 



Using Eq. (A.5) we can calculate the following projec- 
tions of the shear tensor: 

a^n" = En r , a^n" = hi (3n 2 r - l) . (A.8) 

The null geodesic equation can be written 

„ „ dv» 



dX 



0, 



(A.9) 



where A is an affine parameter along the geodesic. Our 
goal will be to evaluate the various components of 
the geodesic equation. For the time component, using 
Eqs. (A.3) and (A.8) we find 



dK ^=-g = ^ + ^(3^-l). (A.10) 



dX 



To extract the remaining components, we can first ex- 
pand the geodesic equation using Eqs. (A.6), (A.3), and 
(A. 10) to obtain 

dn^ 1 1 

= - 7 U M _ + _ 7S ( 3n 2 _ ^ K + 

(A.ll) 

Then we can write the radial component as 



d(v^r^) d( 7 n^r p ) 



dX 



dX 

g?7 dn p 
dX 



(A.12) 



(IT 

n r + 1 —r„ + 1 n^. (A.13) 
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The third term on the right-hand side of Eq. (A. 13) can 
be evaluated as follows: 



, dr 



dX 



(u v + n v ) (A.14) 

7 n e 0% ; „n e ^ (A.15) 

jnle^D^ (A.16) 

\in 2 e 9 3 . (A.17) 



Here expression (A.15) follows from the normalization of 
r M and 7i M , spherical symmetry, and Eq. (A. 3). Expres- 
sion (A.17) uses the normalization of and the defini- 
tion 



(A.18) 



for the "sheet expansion" 9 S , which measures the spatial 
rate of expansion of the spatial congruence r M on each 
comoving time slice. Combining Eqs. (A. 10), (A. 11), and 
(A.17), we finally obtain the radial component of the 
geodesic equation: 



d{v^i 



dX 



= -7 n r H« 



1 2 

2 ' 



(A.19) 



We can obtain the final # M component of the geodesic 
equation in a completely analogous manner to the ra- 
dial component. However, we can do this more easily by 
taking the derivative with respect to A of the identity 
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2 2, 22 

7 n r + 7 n e 



(«%) 2 + (v^Orf 



(A.20) 



After substituting Eqs. (A. 10) and (A.19), the result is 



d{v^ 



^ 2 9 s n r n e . 



(A.21) 



We now have all components of the null geodesic equa- 
tion in covariant 1 + 1 + 2 notation. Our final step will be 
to express these equations in terms of coordinates along 
a null ray. The most natural choice of coordinates is the 
comoving, proper time coordinates implied by the LTB 
metric, Eq. (1). The single nontrivial angular coordinate 
will be taken to be the standard spherical coordinate 9 



{<j> will be constant). The link between the covariant 
geodesic equation and the coordinates x^ is through the 
relation i+ = dx^ j dX. In component form, we have 



vu » = ~dX 



dt 
IX' 
V 



dr 



y/T=KdX 
A9 



Y- 



dX' 



(A.22) 
(A.23) 
(A.24) 



Inserting Eqs. (A.22) to (A.24) into the geodesic equation 
components, Eqs. (A. 10), (A.19), and (A.21), we finally 
obtain the coupled set of ordinary differential equations 



dt 

dX 
dr 

Tix 

d6 
dX 
dH 
dX 2 

dX 2 





dn r 
~dX 



7: 



Y> 



±7Vi 



nt 



Y 



= -7 



Y 7 
d9 

(7a 



K> 



2(1 - K) 

2 F(l - K) 
Y' : 



Y- 



,d9 
dX 



d_ 

dX 

57("?-l) 



3Sn r 



Y 



(A.25) 
(A.26) 

(A.27) 
(A.28) 

(A.29) 
(A.30) 
(A.31) 



Here the expression 9 a = 2\/l — K/Y has been used. 
Equation (A.30) represents the conservation of angular 
momentum. It is straightforward to verify that these 
equations are equivalent to those derived by more famil- 
iar techniques (e.g., [50]). The set (A.25) to (A.31) is 
overdetermined, and experimentation revealed the sub- 
set (16) to (19) (specialized to the decaying mode case 
of K (r) = 0) to be the most robust and numerically effi- 
cient. 
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